GDP = [4.321 4.2352 4.2858 4.3734 4.268 4.2536 4.1162 4.276 4.3617 4.9984];
job = [0.79 0.79 0.81 0.82 0.83 0.83 0.83 0.83 0.84 0.85];
a = polyfit(job,GDP,1);
m = 0.78:0.05:0.85;
y = a(1)*m + a(2);
plot(job,GDP,'.',m,y,'g','LineWidth',2);
axis([0.78 0.85 4 5]);
xlabel('\bfx_{2}');
ylabel('\bfy_{1}');
title('\bfy_{1} = 5.56x_{2} - 0.22')